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Abstract 

A non homogeneous spin chain in the representations {3} and {3*} of A2 is analyzed. 
We find that the naive nested Bethe ansatz is not applicable to this case. A method inspired 
in the nested Bethe ansatz, that can be applied to more general cases, is developed for 
that chain. The solution for the eigenvalues of the trace of the monodromy matrix is given 
as two coupled Bethe equations different from that for a homogeneous chain. A conjecture 
about the form of the solutions for more general chains is presented 



Integrable magnetic chains are interesting physical systems with a rich mathematical 
structure. Homogeneous integrable spin chains with spin 1/2 and higher have been found 
and solved Other interesting systems are non-homogenous chains combining two 

different kinds of spin states in the sites. De Vega and Woynarovich || pioneered the work 
in these systems, whose thermodynamic limit was obtained by de Vega, L. Mezincescu and 
R.I. Nepomechie [[F) (see too H8|-||10||). These authors considered an alternating chain with 
spin 1/2 and spin 1 associated to a su(2) algebra. Following the quantum inverse scattering 
method (QISM) and using the Yang Baxter equation (YBE) |nj, they found integrable 
hamiltonians that present terms coupling pairs of neighboring sites and others coupling 
three neighboring spin sites. The method introduces a monodromy matrix on an auxiliary 
space whose elements are tensors on the state space of the chain, tensorial product of 
the site spaces. The spectrum and eigenstates of the trace of the monodromy matrix are 
obtained by using the Bethe ansatz (BA) method and solving the set of equations, known 
as Bethe equations, derived by this method. 

In a homogeneous system with the associated algebra of higher rank ( ||12|| - [p!5[ ), the 
solutions are obtained by applying the BA method successively (d — l)-times, d being the 
dimension of the spin. This method is know as nested Bethe ansatz (NBA) ||16|| . 

In this paper, we are going to solve a non-homogeneous chain with spin associated to 
the algebra sit (3) and states of the sites alternating between the elementary representations 
{3} and {3*}. In order to obtain hamiltonians associated to alternating chains, we have 
made an extension of the method used in ref. || for systems where P and T symmetry are 
not conserved . The diagonalization of our hamiltonians is not possible by the usual NBA 
method. The reason why is that we have not a reference state, which would be taken as 
pseudovacuum annihilated by the elements under the diagonal of the monodromy matrix, 
to build the eigenstates of the diagonal elements of the monodromy matrix whose trace is 
related to the hamiltonian. 

In the homogeneous case, the eigenstates are built by the action of the elements above 
the diagonal of the monodromy matrix on the pseudovacuum state; this fact makes the 
charges of su(3) be related to the number of applied operators, which can be considered 
as a particle number. In the non-homogeneous case this does not happen anymore, since 
there are site states with charges of different sign which belong to different representations 
and therefore there are states of the chain with the same total charge and different number 
of particles. This makes the usual NBA method non applicable. 
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The method we are going to present is inspired in NBA one and in the method given 
for solutions of 0{2N) symmetric theories in ref. |T7[], but it has important differences with 



them. In a first step, we look for a subspace of states of the chain that are eigenstates of 
the element Ti t ± of the monodromy matrix and are annihilate by the elements T^i (j > 1), 
the whole subspace being invariant under application of T^j > 1). On this subspace, 
we build the eigenstates of the trace of T by the application of operators T±j (j > 1) and 
obtain the eigenvalues in terms of the restriction of T to one less dimension. The procedure 
should be repeated to end with a monodromy matrix of dimension one. As we see, our 
method can be considered as an extension of the NBA one. 

In this paper we start by showing the general features of a non-homogeneous chain 
0. Then develop our method in an alternating chain based in the representations {3} 
and {3*} of su(3), but the method has a general application that we will applied in a 
forthcoming paper. 

The YBE is usually written 

R(6 - B') ■ [t{8) ® t(d')] = [t{9') ® t{B)\ ■ R{0 - $'), (1) 



with 



[taA0)U = R b c d a(0) (2) 



and the product (g> being understood, as usual, in the indices c and d which lie in the site 
space and the • product in the space tensorial product of the two auxiliary spaces in which 
lie the indices a and b of t and all indices of R. Both operators are represented in fig. 1(a). 

2 



Fundamental representation space {3} 
Fundamental representation space {3*} 



[ta,b (6)] c , d = 



(a) 



(b) 



R* M (9) 

7, o v ' 



a V5 p 



(c) 



tUpO)] a ,b = 

(d) 



P 



b 



Fig. 1 



There is another solution to (U) with the same .R-matrix and a t-matrix, with the 
site space in other representation, that we call t and is represented in fig. 1(c) with self- 
explanatory notation. The YBE ([I]) for these operators is represented in fig. 2. 





Fig.2 



As said before, the key of the method is to build the monodromy operator as a product 
of t and t matrices, one per site, in the indices of the auxiliary space. The method could 
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be applied to any distribution of the spin representations in the chain. Here, we are 
considering an alternating chain. So, we have for the monodromy operator r( aW ) 

Tif >(*, a) = 41(0)41,(0 + a) • ■•t^WC-Jfi + «). ( 3 ) 
is graphically expressed in fig. 3(a) 




(b) 

Fig. 3 



This operator, due to (0), fulfills the YBE. 

In the same way we have obtained this system, we can considered an analogous sys- 
tem by interchanging the representations {3} and {3*}. The corresponding R*, t* and t 
operators, represented in fig. 1(b) and (d), satisfy 

R*(9 - 6') ■ [t*{B) ® t*{6')) = [t*(8') ® t*(8)] -R*(6- 9'), (4) 

with 

KMkd = R% d (0). (5) 

The corresponding monodromy matrix, T^ alt \ is represented in fig. 3(b). 
In the most general case, the R -matrix has the following properties 
i) It is PT-symmetric 

R^) = R^j{d). (66) 
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ii) It is unitary 

R c a d b (e)R e c J d (-0)=p(e)S a , e S bJ , (7a) 

R%(0)R^:s(-o)=P*(0)s a ,,s p , v . (76) 

iii) The values for the spectral parameter 6 = are R(0) = cqI and R*(0) = CqI 

iv) A matrix M and a constant r\ exist such that 

R c a ' d b (9)M b , e Rf; d (-e - 2r,)M£ oc 5 a J c , h . (8) 

v) The t-matrices verify 

[ia, b m*Ah,l(- e )kc = p(0)Sa,cS a ,j. (9) 

The hamiltonian is obtained by taking the trace of T^ alt ^ in the auxiliary space 

r<«">(0,a) = T#*>(M), (10) 

and doing 

H = ^^ t \e,a)\ e=0 . (11) 

The successive derivatives of r^ 1 ^ give a complete set of commuting operators that make 
the system integrable. 

Unlike in the homogeneous case, the hamiltonian has pieces that involve two and three 
neighboring sites 

2N-1 2N-2 

iodd iodd 

with 

(/4j+ik/3;&,7 = ?a,c(a)k«fo,7(- Q! )k& (13) 

and 

(^• 5 / + l ii+2 )a,/3, C ;6,7,d = [*a ) e(a)]/9 ) «[*e I d(0)] C) /[t« l7 (-a)]/ ) 6 (14) 

where t means dt/dO. 

In the same way, we can obtain another system by interchanging the representations 
{3} and {3*}. The hamiltonian H so obtained commutes with H and both of them can 
be simultaneously diagonalized. 
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The chain we are considering, the simplest one in a higher rank algebra, combines 
alternating the two elementary representations of su(3). In fig. 1, we associate the solid 



line to the {3} representation and the wavy line to the {3*}. The t matrix is (16 

3 

t(6) 



3 ^ 3 

: sinh (-0 + 7) 22 e i,i ® e ti + sinh ( ~°) 22 ei 



■'J, J 



1,3 = 1 



(15) 



+ sinh (7) exp[(i 



TX 

j - - sign (i 



j))0}eij ® e 



a 



where the superindex s means site space and, in the {3} representation, the are the 
matrices (e^)/^ = 8i,i8j,m ■ 

The t operator can be written in matrix form 



t(9) 



I a 








Vo 





6 



d 











b 



c 







c 


6 













a 













b 



d 








d 





b 






0\ 




c 

b 
a) 



(16) 



with 



a(6)= sinh (-61 + 7), 

6(0) = sinh (^0), 
c(0) = sinh (7)e 5 , 
d(0) = sinh (7)e~2". 



In the same way, we can write 



t(0) 



/a00c0000d\ 
060000000 
006000000 
000600000 
dOOOaOOOc 
000006000 
000000600 
000000060 
cOOOdOOOa/ 



(17a) 

(176) 

(17c) 
(17d) 



(18) 
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with 



a(0) = sinh(^+!), (19a) 

6(0) = sinh (^(0 + 7)), (1%) 

(0+7) 

c(0) = - sinh (7)e~ =~ , (19c) 

d(6>) = - sinh (7)e-^ . (19d) 



The monodromy matrix T aU can be written as a matrix in the auxiliary space, 

/ A(6) B 2 {9) B 3 (0) \ 
T alt (9) = C 2 {9) D 2 , 2 (9) D 2 , 3 (9) , (20) 
\C 3 (0) D 3 , 2 (9) D 3 ,s(0)J 

whose elements are operators in the tensor ial product of the site spaces. 
The YBE for T alt , can be written in terms of its components 

B{9) <g> B{9') = R {2) {9 - 9') ■ (B(9') <g> B{9)) = {B{9') <g> B{9)) ■ R {2) {9 - 9'), (21a) 
A{9)B{9') = g(9' - 9)B{9')A{9) - B{9)A{9') ■ r^(9' - 9), (216) 
D{9) ® B{9') = g(9 - 9'){B{v) ® D{9)) ■ R {2 \9 - 9') - B{9) ® (r< 2 >(0 - 9') ■ D(9')), 

(21c) 

where 



R^(9) 



/l 0^ 

d o 

6 c 

Vo l/ 



.(2) 



(0) 



h- 
a+ 



;(2) 



(0) 



a+ 
a_ 



(22) 



and 



</(*) = 



6(61) 



h+(9) = 



Key 



For the site states, we use the notation 



h-{6) = 



d{9) 

WY 



(23) 



«= | 



d 





u 






(24) 



Inspired in the NBA method, we look for an eigenstate of A that serves as a pseu- 
dovacuum. Firstly we group two neighbor sites and call Wi the subspace generated by the 
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vectors | u, s > and | w, d > in the two site space formed with the 2z-th and (2i + l)-th 
sites. Then, we build the subspace 

O = wi <S> W2 <8> • ■ • <S> wn (25) 

in the total space of states of a chain with 2N sites. 

In a nonhomogeneous chain, we have not a state || s > such that Dij || s >oc 6ij \\ s >. 
For this reason, the NBA method can not be used. Our method, instead, starts with a 
state || 1 >G ft verifying 

A(6) || 1 >= [a(e)] Nsi [b(e)] N Z || 1 >, (26a) 

B i ||l>^0, i = 2,3, (266) 

d\\l>=0, i = 2,3, (26c) 

Di,j || l>eJ2, i,j = 2,3, (26d) 

^3 (A^g) being the number of sites in the representation {3} ({3*}). In our case ./V3 = 
iV 3 * = N. 

Following the steps inspired in the NBA, we apply r-times the B operators to || 1 > 
and build the state 

••,//,.) =B il (fi 1 )®...®B ir (n r )X il „. >ir || 1 >, (27) 

Xii,—,i r being a r-tensor that, together with the values of the spectral parameters 
fix, - ■ ■ , fi r , will be determined at the end. 

The action of A(fj,) and Di^fj,) on \l/ is found by pushing them to the right through 
the Bi 5 (njYs using the commutations rules ( [216, cj) . Two types of terms arise when A and 
Dij pass through S's, the wanted and unwanted terms, similar to the obtained in the 
NBA method. The first one comes from the first terms of ( |216, d. ). In this type of terms 
the A or Di^ and the S's keep their original arguments and give a state proportional to 
\E f . The terms coming from the second terms in ( [216, c| ) are called unwanted since they 
contain Bi((i) and so they never give a state proportional to so, they must cancel each 
other out when we sum the trace of T alt . 

The wanted term obtained by application of A is 

r 

[a(n)] N *$(fj)] N S J] g(fij - y)BiM ® • • • ® B ir (fi r )X iu ..., ir \\ 1 >, (28) 



and the k-th unwanted term 

r 

3 = 1 

•••(g) B(// r ) (2) ® fl( A i fc _ 1 )M<*- 1 )X || 1 >, 

(29) 

M being the operator arising by repeated application of (|21a| ) , 

B(m) ® • • • <g> B(/x r ) = B(/x fe+ i) <g> • • • S(/i r ) <g> • • • ® B(/x fc _i)M ( fe_1 \ (30) 

The application of the operators D^j(jj) to the state \I/(/Z) results more complicated. 
The wanted term results 

r 

[Dkj(M)B il (ni)®...®Bi r (ji r )X iu ... tir || l>] Manted = n^ - ^*) 

i=l (31) 

^Sr 1 ^ - *o • • -i?g?: 2 1,i2 ^ - Ma) • ^Sf^ - ma^h,.,*, n 1 >, 

where the R^'s product is taken in the auxiliary space and has the form 

ti>(„ Tf\ - n(2)ar-i,V n(2)ai,« a R (2)j'.*i _ ( tf) P{t l itf)\ (o \ 
*CM,/*)w=%,a r ■••^ 2 ,a 2 -^(^^ (32) 

The action of Dkj with fe ^ j on || 1 > is not zero. This is the main difference with 
the models that can be solved by NBA. Then, we try to diagonalize the matrix product 

F ( ^ = D W .^={^% (33) 



By taking the terms in (^Tj) with k = j and adding them for k = 2 and 3, we obtain the 
wanted term 

r 

II ^ - (m) ® . . . ® 5 !r (// r )T (2) (|i, (1)4 ... ir || 1 >, (34) 

3=1 

where 

T(2) fa p) = trace(F) = A (2) {fi,jl) + £> (2) (ji,jZ). (35) 
In the same form, the fc-th unwanted term results 



Y[ g(^k - Hj) (B(jj)rM(jji - /Ik)) ® B(fi k+1 ) ® • • • 

V 7 (36) 



3=1 



5(// r ) <g> S(/n) <g> S(^_ 1 )M^- 1) r (2) (^, || 1 > . 



The sum of the wanted terms and the cancelation of the unwanted terms give us the 
relations 

r {2) (fi,fl)X || 1 >= A (2 )(//,/2)X || 1 > (37) 

and 

a (2) o**,aO = tofW* 1 ' II g ~P " ■ (38) 

We must now diagonalize (p7|). 

The tensor X^ ... ^, (ij = 2, 3) lies in a space with 2 r dimensions and || 1 >G 0. Then, 
the vector X || 1 > yields in a space f^ 2 -* with 2 r+Ar dimensions. In this space, we take 
the element 

" 1>(2)= (o) ®""®(o) ® '^ >x I^>n, ( 39 ) 

which is annihilated by C*- 2 -* (//, (Note that the operators a, /3, 7 and 5 of (|3~2| ) act on 
first part of || 1 >^ 2 - ) and the operators Dij on the second part ). The application of the 
operators and gives 

A^(^fl) || 1 >( 2 )= [&(/*»)]"" U»k)] N; II 1 > (2) , (40a) 
Z^ 2 W) || 1 >( 2 )= n^^faWl^rW^ II 1 > (2) , (406) 

The important fact is that .F(/U, /7) verifies the YBE with the matrix given in 

(H, 

(jt-y!) /Z) <g> = [iV, # ® F(/x, A?)]^ (2) (// - //) , (41) 

which, in a second step, permits us to solve the system. From this equation, we obtain the 
commutation rules 

A^{y) • B( 2 V) = g(y' - fi)B^(y') • ^ 2) (y) - h+tf - ^)5 (2) M • A< 2 V), (42a) 
D&(y) ■ B^(y') = g(fi - /*')B (3 V) ■ DV>{ji) - h+Q* - n')B^{y) • D^(y'). (426) 



In this second step, we build the vector 

vl>( 2 )(A,^)= J B( 2 )(A 1 ,^)--- J B( 2 )(A s ,^) ||1>( 2 ) 
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(43) 



The action of A^\X, ft) on # (2) gives the wanted term 

s 

[b(X)] N -[b(X)] N S]Jg(X i -X)B( 2 \X 1 ,jl)...B^(X s ,fl) \\1>^\ (44) 



i=l 

and the k-th unwanted term 



(45) 



-h+(X k - X)lb(X k )} N *[b(X k )] N 3 J] g(Xi - X k )B&(X, jl)B^(X k+1 , ft)... 

i = l 

. . .S (2) (A fc _i, /2) ||1>( 2 ). 
In the same form, the action of D^ 2 \X, p) on \I/( 2 ) gives the wanted term 

S T 

[b(X)] N3 [a(X)] N * J\g(X - A,) n 1 B< 2 >(Ai, ft . . .5< a >(A., # II 1 > (2) , (46) 
i=i j=i ^ Afc ^ 

and the k-th. unwanted term 

S T 

-h- (A - X k ) [b(X k )] N * [a(X k )] N * J] g(X k - A,) j[ _J_fl( 2 > (A, p) B^ (A fc+1 , /2) . . . 



...£( 2 )(A fe _ l5 /2) ||1>^. 

(47) 

The cancelation of the unwanted terms and the sum of the wanted terms, give us the 
equations 

a(Afc) ] ^ A 1 _ A ff(A t - A fc ) _ 
6(A fc )J f = 1 1 ^(A fc -/x i )"M^(A & -A i ) fc - i ---' s ^ 

and 

s 

A( 2) (A**, AO = II ^ - Mi) [&(A*)] [S(A**)1 ^ • (49) 

i=l 

Then, by comparing equations fl38| ) and (|4"9| ) and calling g{9) = a(9)/b(6), we obtain 
the coupled Bethe equations 

[g(Xk)] N * = n 9(X k - H) II - \ k \ , (50a) 

s 
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and the eigenvalue of the trace of T^ alt ^ 
m = [a(fJ,)) N nK»)} N °il9(H-»)+ 



3=1 

r 



7=1 9\r ^j, 



3=1 



1=1 1=1 



(51) 

that is the solution to our problem. 

The solution we have obtained, permits us to conjecture the form of the solution for 
a non homogeneous chain combining the two elementary representations of an algebra A n . 
Each elementary representation introduces a function g and g (that we can call source 
functions) . Such solution will have n sets of Bethe equations ( the same number of basic 
representations of A n ) . The first and last sets will have in the first member their respective 
source functions powered to the number of sites of each representation and in the second 
member a product of source functions similar to those in the right hand side of eqns. 
( [50a, Vj) . The left hand sides of the rest of the sets of equations are the unity, since there 
is not source function. For other representations, we expect a similar scheme with source 



functions related to the source functions of the elementary representations . Since our 
method has a wide applicability, we will study these points in a forthcoming paper. 
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